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We introduce two possible ways of defining effective constraints of quantum systems and applied 
this effective constraint method to models of WDW Quantum Cosmology and Loop Quantum Cos- 
mology. We analyze effective Hamiltonian constraint on both second and third order and calculate 
Hubble parameter as well as modified Friedmann equation of each model. Then we compare with a 
special case using coherent state. It shows that this method is reasonable and as before the classical 
' Big Bang singularity is replaced by a quantum bounce in Loop Quantum Cosmology. 
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. I. INTRODUCTION 



In quantum systems, there are many applications which are placed in regions close to classical behavior, where 
quantum properties must be taken into acount. So it is often helpful to work with equations of classical type, i.e. 
sets of ordinary differential equations for mechanical systems, which are amended by correction terms resulting from 
quantum theory. Then from a mathematical point of view, the question becomes how well the behavior of a wave 
function subject to a partial differential equation can be approximated by finitely many variables subject to a system 
of coupled but ordinary differential equations. 

One useful example is the effective equations of motion of quantum mechanics. Quantum mechanics can be for- 
mulated into an intrinsically geometric fashion[l]. Based on this formulation, we can develope a method of effective 
equations of motion which is a trustworthy tool to avoid some of the mathematical and conceptual difficulties of quan- 
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of outright quantum states, and they can often be derived in a manageable way [2]. Effective equations are not only 



^ . quantum corrected classical equations, but also provide direct solutions for quantum properties such as expectations 



values or fluctuations. While semiclassical regimes play important roles in providing useful approximation schemes, 
effective equations present a much more general method. 

One of the key issues in quantizations of fundamental theories is the determination of physical observables. They 
must satisfy the constraint equations and be invariant under gauge transformations. For canonical quantum theories, 
solving constraints is traditionally done at the state level: one constructs a physical Hilbert space of states annihilated 
by the constraint operators and equipped with an invariant inner product. Explicit constructions can be done in some 
[CI. ■ special cases by different methods. Since explicit derivations are possible only in specific cases, it is not always 
clear whether the results are generic or mere artefacts of the simple models used. It is therefore important to have 
approximate methods for a wider range of cases, or at least to be able to perturb around known solvable ones while still 
' ensuring that the constraints are solved and the observables are gauge invariant. It turns out that such perturbation 
, schemes are most feasible if one deals with the observables directly, such as expectation values. This procedure gives 
fNj ■ rise to canonical effective equations and constraints [2, 4]. 

t-H | So for constrained systems, we can develop an effective constraint formalism parallel to that of unconstrained 
systems. Its advantages are that (i) it avoids directly writing an integral (or other) form of a physical inner product, 
which is instead implemented by reality conditions for the physical variables; and (ii) it directly provides physical 
quantities such as expectation values and fluctuations as relational functions of internal time, rather than computing 
a whole wave function first and then performing integrations. These advantages avoid conceptual problems and some 
technical difficulties in solving quantum equations. They can also bring out general properties more clearly, especially 
in quantum cosmology [11, 12]. 

One of the most fundamental issues in modern physics is to establish the theory of quantum gravity. Among 
various approaches, the nonperturbative quantization is a very important candidate. Based on the ADM Hamiltonian 
formalism of general rclativity(GR) and the Dirac's generalized Hamiltonian quantization method, Wheeler and 
DeWitt proposed a wave functional description of gravitational field. In this quantum geometrodynamical approach, 
the whole universe is described by a wave function which is defined on superspace, the space of all 3-metrics and matter 
field configurations. The dynamics is encoded in the second-order differential Wheeler-DeWitt(WDW) equation. 
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There is an alternative background independent approach which is the loop quantum gravity(LQG) [15-18]. The 
starting point of LQG is the Hamiltonian connection dynamics of GR rather than the ADM formalism. In this 
framework, GR looks like a gauge held theory with SU(2) as its internal gauge group. By taking the holonomy of 
su(2)-connection A l a and flux of densitized triad E b as basic variables, the quantum kinematical framework of LQG has 
been rigorously constructed, and the Hamiltonian constraint operator can also be well defined to represent quantum 
dynamics. Moreover, a few physically significant results, especially the resolution of big bang singularity, have been 
obtained in the minisuperspace models of loop quantum cosmology (LQC) [19-23]. The quantum bounce replacement 
of big bang and its properties are being studied from different prescriptions of LQC [24]. Effective equations were also 
derived in isotropic models [25, 26], which predict evolution of universe with quantum corrections and shed new lights 
on the singularity resolution. 

In both WDW quantum cosmology and LQC, by coupling with a massless scalar field, the Hamiltonian constraint 
equations can be reformulated as Klcin-Gordon-like equations, where the corresponding gravitational Hamiltonian 
operators, as multiplications of several self-adjoint operators, are non-symmetric in the kinematical Hilbert space. 
The above models of quantum cosmology provide a good arena to launch the desired investigation. Hence we will 
study the effecive constraints in spatially flat isotropic FRW cosmology coupled with a massless scalar field in both 
WDW and loop quantum cosmology. 

We give a brief introduction of effective equations of motion for unconstrained system in section II and two definitions 
of effective constraints in section III. Then in the major part of this paper, section IV, we analyze the second and 
third order effective constraints in WDW Quantum Cosmology and Loop Quantum Cosmology. Finally, we compare 
the result with the one in a special case using coherent state in the discussion section. 

II. GENERAL FORMALISM 

We first review the setup of effective equations for unconstrained Hamiltonian systems [2, 3], which we will generalize 
to systems with constraints in the following section. 

We describe a state by its moments rather than a wave function. This has the immediate advantage that the 
description deals directly with quantities of physical interest. In terms of expectation values, fluctuations and all 
higher moments, this structure takes the form of an infinite dimensional phase space whose Poisson relations are 
derived from the basic commutation algebra. Dynamics is determined by a Hamiltonian on this phase space. As a 
function of all the phase space variables it is obtained by taking the expectation value of the Hamiltonian operator 
in a general state and expressing the state dependence as a dependence on all the moments. Thus, the Hamiltonian 
operator determines a function on the infinite dimensional phase space which generates Hamiltonian evolution. 

For an ordinary quantum mechanical system with canonical basic operators q and p satisfying [q,p] = ih, we have a 
phase space coordinatized by the expectation values q :— (q) and p :— (p) as well as infinitely many quantum variables 

G"> b :=((p-(p)r( q -(q)) b ) Weyh (2.1) 

for integer a and b such that a + b > 2, where the totally symmetric ordering is used. For a + b = 2, for instance, 
this provides fluctuations (Aq) 2 = G°< 2 = G qq and (Ap) 2 = G 2 <° = G pp as well as the covariance G M = G qp . As 
indicated, for moments of lower orders it is often helpful to list the variables appearing as operators directly. The 
symplectic structure is determined through Poisson brackets which follow by the basic rule {A, B} = — ih~ 1 ([A, B]) 
for any two operators A and B which define phase space functions A := (A) and B := (B). Moreover, for products of 
expectation values in the quantum variables one simply uses the Leibniz rule to reduce all brackets to the elementary 
ones. General Poisson brackets between the quantum variables then satisfy the formula 

{G a ' b , G°' d } = -£ (-^y 

r s=0 j k \J J V / V / \J J 

+adG a - 1 > h G c > d - 1 - bcG^^G"- 1 ^, (2.2) 

where the summation of j and k is over the ranges < j < min(a, d) and < k < min(6, c), respectively. For low 
order moments, it is easier to use direct calculations of Poisson brackets via expectation values of commutators. For 
second order moments, it is easy to get 

{G^G^j^G ' 2 , {G^G^j^G 1 ' 1 , {G 1 ' 1 , G 2 < } = 2G 2 <°. (2.3) 

This defines the kinematics of the quantum system formulated in terms of moments. 
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Dynamics is denned by a quantum Hamiltonian derived from the Hamiltonian operator by taking expectation 
values. This results in a function of expectation values and moments through the state used for the expectation value. 
By Taylor expansion, we have 



H Q (q,p,G a > b ) = (H(q,p) Woyl ) = (H{q + (q-q),p+(p-p)) Wcyl ) 

1 d a+b H(q 
a!6! dp a dq 



a=0 6=0 

where we understand G a,b = 0ifa + &<2 and H(q,p) is the classical Hamiltonian evaluated in expectation values. 
As written explicitly, we assume the Hamiltonian to be Weyl ordered. If another ordering is desired, it can be reduced 
to Weyl ordering by adding re-ordering terms. 

Having a Hamiltonian and Poisson relations of all the quantum variables, one can compute Hamiltonian equations 
of motion q = {q : Hq} : p = {p,Hq} and G a ' b = {G a ' b ,HQ}. This results in infinitely many equations of motion 
which, in general, are all coupled to each other. This set of infinitely many ordinary differential equations is fully 
equivalent to the partial differential equation for a wave function given by the Schrodinger equation. In general, 
one can expect a partial differential equation to be solved more easily than infinitely many coupled ordinary ones. 
Exceptions are solvable systems such as the harmonic oscillator or the spatially flat quantum cosmology of a free, 
massless scalar field [7] where equations of motion for expectation values and higher moments decouple. In quantum 
cosmology, this is developed in [9, 10, 12]. Moreover, semiclassical and some other regimes allow one to decouple and 
truncate the equations consistently, resulting in a finite set of ordinary differential equations. This is easier to solve 
and can be exploited to avoid conceptual problems especially in the context of constrained systems. 



III. EFFECTIVE CONSTRAINTS 



For a constrained system, the definition of phase space variables (2.1) in addition to expectation values of basic 
operators is the same. For several basic variables, copies of independent moments as well as cross-correlations between 
different canonical pairs need to be taken into account. A useful notation, especially for two canonical pairs (q, p; qi,pi) 
as we will use it later, is 

Gci = ((P-P) a (<l <?)Vi -Pi) c (4i <?i)Vc yl . (3.1) 

Also here we will, for the sake of clarity, sometimes use a direct listing of operators, as in G qq — Gq'q = (Aq) 2 or the 
covariance G q 1 — G^'J, for low order moments. 

We assume that we have a single constraint G in the quantum system and no true Hamiltonian; cases of several 
constraints or constrained systems with a true Hamiltonian can be analyzed analogously. We clearly must impose the 
principal quantum constraint CQ(q,p,G a ' b ) := (C) = since any physical state |\&), whose expectation values and 
moments we are computing, must be annihilated by our constraint, C\^} — 0. Just as the quantum Hamiltonian Hq 
before, the quantum constraint can be written as a function of expectation values and quantum variables by Taylor 
expansion as in (2.4). However, this one condition for the phase space variables is much weaker than imposing a Dirac 
constraint on states, C\^) = 0. In fact, a simple counting of degrees of freedom shows that additional constraints 
must be imposed: One classical constraint such as C — removes a pair of canonical variables by restricting to the 
constraint surface and factoring out the flow generated by the constraint. For a quantum system, on the other hand, 
we need to eliminate infinitely many variables such as a canonical pair (q,p) together with all the quantum variables 
it defines. Imposing only Cq = would remove a canonical pair but leave all its quantum variables unrestricted. 
These additional variables are to be removed by infinitely many further constraints. 

Here we introduce two ways of defining quantum constraints. The first way is to define constraints with operators 
of non-symmetric ordering. There are obvious candidates for these constraints. If G|^) = for any physical state, 
we can define infinitely many quantum constraints 

C {n) ._ ^ = Qj (3.2) 
C flL) ■= </(<?'P)G">=0, (3.3) 

for positive integer n and arbitrary phase space functions ,f(q,p). All these expectation values vanish for physical 
states, and in general differ from each other on the quantum phase space. For arbitrary f(q,p), there is an uncountable 
number of constraints which should be restricted suitably such that a closed system of constraints results which 
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provides a complete reduction of the quantum phase space. The form of functions f(q,p) to be included in the 
quantum constraint system depends on the form of the classical constraint and its basic algebra. 

We thus have indeed infinitely many constraints, which constitute the basis for our effective constraints framework. 
This is to be solved as a classical constrained system, but as an infinite one on an infinite dimensional phase space. 
An effective treatment then requires approximations whose explicit form depends on the specific constraints. At this 
point, some caution is required: approximations typically entail disregarding quantum variables beyond a certain 
order to make the system finite. Doing so for an order of moments larger than two results in a Poisson structure 
which is not symplectic because only the expectation values form a symplectic submanifold of the full quantum phase 
space, but no set of moments to a certain order does. We are then dealing with constrained systems on Poisson 
manifolds such that the usual countings of degrees of freedom do not apply. For instance, it is not guaranteed that 
each constraint generates an independent flow even if it weakly commutes with all other constraints which would 
usually make it first class. Properties of constrained systems in the more general setting of Poisson manifolds which 
need not be symplectic are discussed, e.g., in [13]. 

Some of quantum constraints in (3.3) are defined as expectation values of non-symmetric operators, thus implying 
complex valued constraint functions. We specifically do not order symmetrically in (3.3) because this would give rise 
to terms where some q or p appear to the right while others remain to the left. This would not vanish for physical 
states and therefore not correspond to a constraint. This may appear problematic, but one should note that this 
reality statement is dependent on the kinematical inner product used before the constraints arc imposed. This inner 
product in general differs from the physical one if zero is in the continuous part of the spectrum of the constraint 
and thus reality in the kinematical inner product is not physically relevant. Moreover, in gravitational theories it is 
common or even required to work with constraint operators which are not self-adjoint [14], and thus complex valued 
constraints have to be expected in general. For physical statements, which are derived after the constraints have been 
implemented, only the final reality conditions of the physical inner product are relevant. 

This physical reality can be implemented effectively: We solve the constraints on the quantum phase space, and 
then impose the condition that the reduced quantum phase space be real. We will see explicitly that complex- valued 
quantum variables on the unconstrained phase space are helpful to ensure consistency. In parallel to Hilbert space 
notation, we call quantum variables (2.1) on the original quantum phase space kinematical quantum variables, and 
those on the reduced quantum phase space physical quantum variables. Kinematical quantum variables are allowed 
to take complex values because their reality would only refer to the inner product used on the kinematical Hilbert 
space. For physical quantum variables in the physical Hilbert space as usually defined, on the other hand, reality 
conditions must be imposed. 

It may seem obvious how to avoid the question of reality of the constraints altogether by using quantum constraints 

defined as G c "^^' p ' = (C n f(p,q))weyi such as G" 7 " 9 and G c " p with the symmetric ordering used as in (2.1). This 
is second way of defining the effective constraints which are symmetric. Here, the symmetric ordering contained in 
the definition of quantum variables must leave C intact as a possibly composite operator, i.e. we have for instance 
G Cp — \ {Cp + pC) — Cp independently of the functional form of C in terms of q and p. Otherwise it would not be 
guaranteed that the expectation value vanishes on physical states. We could not include variables with higher powers 
of q and p, such as G c pp as constraints because there would be terms in the totally symmetric ordering (such as 
pC n p) not annihilating a physical state. But, e.g., G Cp understood as \{Cp 2 + p 2 C) — Cp 2 would be allowed. The 
use of such symmetrically ordered variables would imply real quantum constraints. 

However, this definition of quantum constraints may seem problematic: The constraints would not form a closed 
set and not even be first class. We have, for instance, 

{ G c"/( 9 ,p) ;G c m ff ( g ,p)| = A_^c n f + fC n ,C m g + gC m ]) 

-^([C^f + fC^C™]) - —([C n f + fC n ,g]) 

f r<n 
-^ Cm > &n 9 + aC m ]) - 2^<[/, C m g + gC m \) + {C n f, C m g}. 

The first commutator contains several terms which vanish when the expectation value is taken in a physical state, but 
also the two contributions [C n ,g]C m f and fC m [C n ,g] whose expectation value in a physical state vanishes only if / 
or g commute with C . This would require quantum observables to be known and used in the quantum constraints, 
which in general would be too restrictive and difficult. In some models, we may lucky to find out enough quantum 
observables. 

By contrast, the quantum constraints defined in the first way do form a first class system: We have 



[fC n , gC m ] = [/, g]C n+m + f[C n ,g]C m + g[f, C m ]G 



(3.4) 
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whose expectation value in any physical state vanishes. Thus, using these constraints implies that their quantum 
Poisson brackets vanish on the constraint surface, providing a weakly commuting set: 

{C (p) iC (rn) } = L^ffr » 0. (3.5) 

Constraints thus result for all phase space functions f(q,p), but not all constraints in this uncountable set can be 
independent. For practical purposes, one would like to keep the number of allowed functions to a minimum while 
keeping the system complete. Then, however, the set of quantum constraints is not guaranteed to be closed for 
any restricted choice of phase space functions in their definition. If cj™^ and Cg™* are quantum constraints, closure 

requires the presence of (for n > 2), cj"p m g j and C g ^ m ^ as additional constraints according to (3.4). This 

allows the specification of a construction procedure for a closed set of quantum constraints. As we will see in examples 
later, for a system in canonical variables (q,p) it is necessary to include at least Cq 1 ^ and Cp m ^ in the set of constraints 
for a complete reduction. With = ihC^ n \ the first new constraints resulting from a closed constraint algebra 

add nothing new. However, in general the new constraints C^ m p j and C^ m ^ will be independent and have to 
be included. Iteration of the procedure generates further constraints in a process which may or may not stop after 
finitely many steps depending on the form of the classical constraint. 

Although many independent constraints have to be considered for a complete system, most of them will involve 
quantum variables of a high degree. To a given order in the moments it is thus sufficient to consider only a finite 
number of constraints which can be determined and analyzed systematically. 



IV. EFFECTIVE CONSTRAINTS OF QUANTUM COSMOLOGY 

The Hilbert-Einstein action for these models are given by 

S = ^ J d^x^—gR + ± J .l'.r x — ; r\: Kll ,: K! ,. (4.1) 

In this spatially flat model, we fix a space-like sub-manifold S, which is topologically R 3 and equipped with Cartesian 
coordinates x % (i = 1, 2, 3), and a fiducial flat metric °q a b given by °q a bdx a dx h = dx\+dx\ +dx\. The physical 3-metric 
q a b is then determined by a scale factor a satisfying q a b = a 2 °q a b- It is convenient to introduce an elementary cell V 
and restrict all integrations to this cell since the spatial slice is non-compact. The volume of V with respect to °q a b 
is denoted as V and the physical volume is V = a 3 V . Then the geometrical pair (a,p a ) can be used as canonical 
variables, where the conjugate momentum satisfies p a = —j^kaa. 

In order to study the WDW quantum cosmology and LQC on the same footing, we employ the new canonical 
variables (A a ,Ef) in both theories. Due to the homogeneity and isotropy, we can fix a set of orthonormal cotriad and 
triad (°oj a ° e") compatible with °q ab and adapted to V. Then the cotriad uj a which are orthonormal with respect to 
physical metric q a b can be written as uj a — xa°u l a , where \ = 1 if w i h as the same orientation as the fiducial °uj a and 
X = — 1 if the orientation is opposite. The basic canonical variables take the simple form 

A\ = cVo k X, E? = p^qVo 1 °e». (4.2) 

The dynamical variables are thus reduced to (c, p) with the Poisson bracket: {c,p} = 8irG"f /3, where 7 is the Barbero- 
Immirzi parameter. Following the /2-scheme of "improved dynamics" [21], the regulator p, used in holonomies is given 
by p, = y/A/\p\, where A = 2y/3iT"f£ 2 P is a minimum nonzero eigenvalue of the area operator. In order to do the 
semiclassical analysis, it is convenient to introduce new dimensionless conjugate variables: 

c 2nG-/ p a -= sgn(p)|p|i = eV> 3 

■ 2|p|i 3V a a 2 ' 2ttG 7 2ttG 7 ' 1 ' ' 

with the Poisson bracket b, v = 1, where the Planck length £p is given by £ 2 P = Gh. From the matter part of action 
(4.1), we can get the momentum of <j) as p^ — p*<j) = a 3 V <j> and the Poisson bracket: {0,_p</,} = 1. 
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A. Effective Constraints of WDW Quantum Cosmology 

Now we apply the effective constraint method to WDW quantum cosmology considering k = FRW cosmology 
model. We start from the classical Hamiltonian constraint which can be obtained from the action (4.1), 

After quantization, we have the corresponding constraint operator 

C = --bub + — ^ — = Aviv- 1 - Bbub, (4.5) 
7 47TG-7 v v 

where A = 4 J~ Gl , B = ^, and the operator ordering of the first term is to ensure the constraint operator's Hermiticity. 
As usual, basic operators satisfy commutators [<j>,p<t,\ = ih and [b, v] = ih. As in this method, moments are defined as 

Gci = <(* ") a ( b b )\h - p^) c (4> <t>) d )w ey i- (4.6) 

A special class of semiclassical states, such as Guassion states, allows these moments to have the hierarchy 
G°'J = 0{h^ a+b+c+d ^ 2 ). Based on this hierarchy, we can implement a semiclassical expansion, and for semiclas- 
sical approximation, we keep only the moments to a certain order. From now on, we suppose there exists such a 
semiclassical state, on which we calculate all the expectation values. 

We first consider the non-symmetric constraints cj"] p ) := (f(Q,p)C n ). In this model, a complete and first class 
set of quantum constraints can be choosed as 

C Q = (C), C v = {{v - v)C) = {vC) - vC Q , C b = ((b-b)C) = (bC)-bC Q , 

c P , = m-p*)C!) = {p^C) - P4> C Q , c = ((4> - 4)6) = (4>C) - 4C Q . (4.7) 

Some of the constraints are not Weyl ordering, but they only differ by some terms of order of h. 

(v 2 b 2 ) w = vbvb + ih{vb) w + 0(H 2 ), {v-^w = hi)- 1 + \ih v - 2 + 0(h 2 ), 

(vb 3 ) w = b 2 vb - \itib 2 + 0(h 2 ), (p 2 4>)w = # 2 - ihp^. (4.8) 

We first re-order the constraints, then Taylor expland them into functions of expectation values and moments. Con- 
sidering properties of semiclassical state, for second order we only keep terms to the order of H, we thus have 

C Q = Aplv- 1 - Bvb 2 + Av^G^ - 2Ap^- 2 G\'° + Ap 2 v- 3 G 2 Q ° - 2BbG^ - BvG° % (4.9) 
C v = 2Ap 4 ,u- 1 G\% - {Ap 2 v- 2 + Bb 2 )G 2 <° - 2BvbG 1 [ 1 + ihBvb, (4.10) 

C b = 2Ap 4 ,u- 1 G° 1 $ - {Aplv- 2 + Bb 2 )G^ - 2BvbG% - (UhAp 2 ^- 2 + UhBb 2 ), (4.11) 

C p , = 2Ap 4> v- 1 G 2 -(Ap 2 v- 2 +Bb 2 )G 1 1 -2BvbG 1 % (4.12) 
Q = 2Ap^v- 1 G {° 1 -{Ap 2 v- 2 + Bb 2 )G 1 o 1 -2BvbG \ 1 1 +ihAp l pv- 1 . (4.13) 

Here, we can see that is decoupled from other constraints. 

For consistency, the set of all effective constraints should be first class if we start with a single classical constraint 
or a first class set of several constraints. Now we check all Poisson brackets between every two effective constraints. 

{C Q ,C V } = -2BbC v - 2BvC b « 0, {C Q , C } = -2Av~ x C^ + 2A V4> v- 2 C v « 0, 
{C Q , C b } = 2A V4> v- 2 C V4i - 2Bv- x b 2 C v + 2BbC b « 0, {C Ql C P J = 0, 
C b } = -2Ap^v- x C v ^ + 4Bb 2 C v + 4Bvb « 0, {C v , C P J = 2BvbC Pit « 0, 
{a, = -2Ap^C v + 2BvbC 4> w 0, {C 6 , Cj,} = -2A V<i> v- x C b - 2Bb 2 C 4 , w 0, 
{C bl C P J = -2Bb 2 C p ^ w 0, {C p ,, C } = -AAp^C^ + 2Bb 2 C v + 2BvbC b « 0. (4.14) 

It shows that this effective constraint system is indeed a first class system. 

At this point, a further choice arises: we need to determine which variables we want to solve in terms of others 
which are to be kept free. This is related to the choice of time in a deparametrization procedure. Here, we view <j> 
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as the time variable which is demoted from a physical variable to the status of an evolution parameter. Notice that 
time is chosen after quantization when dealing with effective constraints. 

Having made a choice of time, a complete deparametrization requires that all quantum variables of the form G^'^ 

with a 2 7^ or b 2 7^ be completely constrained or removed by gauge. Only quantum variables Gq q are allowed 
to remain free, and must do so without any further restrictions. To second order, the deparametrized system has 
2 + 3 = 5 variables; the parametrized theory has 4 + 10 = 14. We begin by eliminating quantum variables in favor 
of the variables associated with the canonical pair (<fi,p$) only. From the fact that, on the one hand, G^'q, G^'J and 

G ' 2 should satisfy the uncertainty relations and thus cannot all vanish but, on the other hand, are not present in the 
unconstrained system, we expect at least one of them to be removed by gauge. 

P4> = A-iBivb-^B-^-'b-'G^ + v-'Gl* 

A A -h B h u -i bG lfi + A-iBiG^l + l -A^B^b- l Gl% (4.15) 
^B^Gol-^ihA-^BK, (4.16) 

?BivG° ' 2 + ^ihA-iBh, (4.17) 
^vbGl'l + A^Bv^G^, (4.18) 
iBivG° >{-^ih. (4.19) 

The last one is not yet completely expressed in terms of moments only of (q,p). The remaining moments of (0,p^,) 
are not constrained at all, and thus must be eliminated by gauge transformations. To summarize, three expectation 
values are left unconstrained, one of which should be unphysical; six second-order variables are unconstrained, three 
of which should be unphysical. Notice that there is no contradiction to the fact that we have four weakly commuting 
(and independent) constraints but expect only three variables to be removed by gauge. These are constraints on 
the space of second order moments, which, in this truncation, as noted before do not have a non-degenerate Poisson 
bracket (although the space of all moments has a non-degenerate symplectic structure). Weak commutation then 
does not imply first class nature in the traditional sense (see e.g. [13]), and four weakly commuting constraints may 
declare less than four variables as gauge. While the constraints as functionals are independent, their gauge flows may 
be linearly dependent. 

We also emphasize that gauge flows generated by quantum constraints on the quantum phase space play important 
roles. To explicitly account for the unphysical degrees of freedom, we consider the gauge transformations generated 
by the constraints. In our method, Cq is treated as the Hamiltonian constraint which generates the internal time 
evolution. Thus in second order, we only consider gauge flows generated by the four constraints C p<p , G„, Cb and C^ 
but not Cq. After analyzing the gauge transformation, we finally get two observables 

Ox = ub + Gl'l - ^ih, 

0-2 = b 2 G 2 Q '° + 2vbGl'l + » 2G °o 2 o- ( 4 - 2 °) 
And any gauge fixing should not affect these observables. Now we can choose a gauge fixing 

G^ = o, G ;? = G°;l = o, =>G°;° = -±ift, (4.21) 

and let Gq'q, Gq J and Gq'q to be free as physical variables as well as 4> to be free as internal time parameter. 
After substituting the variables required, we can obtain the effective Hamilton constraint, 

C Q = A v \v- X - Bvb 2 - 2BbG 1 ' 1 + ihBb. (4.22) 

It worth noting that this effective Hamilton constraint Cq contains an imaginary term. At the first sight, it seems 
problematic. But when we think it twice, this manifests the non-commutative property of quantum operators. It 
can be expressed into a Moyal *-product inspired by non-commutative geometry. We can see from the derivation of 
equation (4.22) that because of using equation (4.10), it results in an expectation value of a non-symmetric operator, 
and then leads an imaginary term. In the classical constraint, it contains a product of v and b 2 , but after quantization, 



G{'° = A-$B*bGl$+A- 

G% = A-*B*bG$+A- 
G°;° - A- 1 Bb 2 G 2 ° + 2A- 
G?;? = A-iBhG^ + A- 
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they are replaced by the corresponding operators. Considering the non-commutative property of operators, we can 
change the usual product into a Moyal *-product, 

v * b 2 = ue f(Kdt-Kdt) b 2 = vh 2 _ m ( 4 23 ) 

Here we can see the imaginary term is coincident with the imaginary term in Cq . So Cq can be expressed as 

C Q = A v \v- Y -Bv*b 2 - 2BbG 1 Q j ) . (4.24) 

This method is correspoding to a procedure of deparameterization of <f>, 

V \ = A'^BiyH 2 + 2vbG 1 ' 1 - ihvb). (4.25) 

Then we can obtain a true Hamiltonian p^ = ±H , 

H = A-$Bt(yb + G$-*)=A-iBi(v*b + G$). (4.26) 
The dynamics of the system can also be determined by a struture of Moyal *-product, 

f(b,v):=±;(f*C Q -C Q *f). (4.27) 

Since we finally get the effective Hamiltonian constraint Cq , we can use it to generate the evolution of p^ and <j>, 

P* = {P*. C Q } = 0, = {</>, C Q } = 2Au~%. (4.28) 
As expected, (f> can be treeted as internal time generating time evolution. So the time evolution of v is, 

v = {u, C Q } = ^{u *Cq-Cq*v) = 2Bvb + 2BGl:l ~ ihB - ( 4 - 29 ) 



Then we can calculate Hubble parameter, 



4B 



2 



H o = ^ = — ^ + 2 ^ lbG l:l - ■ (4.30) 

In order to eliminate b, we use the effective Hamiltonian constraint again. From Cq = 0, we can have 

b 2 = AB- x V \v- 2 - 2v- 1 bGl A {> + ihv- x b. (4.31) 
Put it into Hubble parameter (4.30), keeping terms to order of h, we yeild the effective Friedmann equation, 

H* =^[AB-ty->]~p. (4-32) 

where p — 2 (2-nG~fv) 2 1S tne ener gy density. This is the same as classical cosmology. There exist a Big Bang singularity. 

In the above treatment, the non-symmetric operators in definition of constraints may give rise to imaginary terms 
in effective constraints. To avoid this imaginary terms, we may change to choose symmetric constraints. But one 
crucial problem is that for this system to be close, we should find out observables before implementing calculation. 
In this system, our goal is to obtain the principle constraint Cq only with physical variables. We may not be able to 
find out enough observables to define constraint at each order and then to solve for all non-physical variables. But 
only for the goal we just indicated, we only need to solve some of the non-physical variables. 

The principle constraint is still 

C Q = Ap\v- X - Bvb 2 + Av^G^ - 2Ap^- 2 G{' a Q + Ap 2 v- 3 G 2 Q fi - 2BbG^ - BisG° ' 2 . (4.33) 
From the constraint operator (4.5), we can find that there is only one observable p$. Then we can define 

C P , = \(p*C + Cp4>) - P4> C Q = 2A P ^- 1 G° 2 ° - {Ap\v- 2 + Bb 2 )G\l - 2BvbG%. (4.34) 
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Obviously one constraint is not enough, and we need to define constraints. We analyze the definition of these symmetric 
constraints in detail. And finally we can find that for this constraint system to be close, we can only include one 
constraint defined from a non-observable. And we choose 

C v = \{i>C+ Cv) - vC Q = 2Ap 4> v- 1 G{'° - (Ap\v~ 2 + Bb 2 )G 2 ^ - 2BvbG 1 j ) . (4.35) 

Now we check the closure of this system. 

{c„,c p j = {^{od + dty-vCQ^fod + dp^-p+CQ} 

= ^{[vC + Cv^C + Cp*}) = ±{ud[d,P4] + [d,P4,]dt>) = 0. (4.36) 

Also we can have a weak constraint, bv and vb are weakly commuted with constraint operator, 

[bv, C] = [vb, C] = -itlC. (4.37) 

Then we can define another constraint, 

Cbu = ^{bv + vb)C+^C{bv + vb))-bvC Q -bC v 

= 2Ap^ 1 -2Bvb 2 G 1 [ 1 -2Bv 2 bG^ 2 . (4.38) 
Now we check whether the Poisson bracket between C v and Cbv is zero, 

{C v ,C bv } = {^{vC + Cv)-vC Q ,^{bvC + vbC + Cbv + Cvb)-bvC Q -bC 1J } 

= ^7z([i>C + Cv,bvC + vbC + Cbv + Cvb]) 

= ^{[C, bv]Cv + [C, vb]Cv + vC[C, bv] + vC[C, vb]) 

= — (iHC 2 v + iHC 2 v + ifivC 2 + ifivC 2 ) = 0. (4.39) 
oin, 

Then we can solve for non-physical observables, 

P(j> = A-iBivb-^B-^-'b-'G^ + v-'G 1 ^ 

~A-iBiv-HG 2 fi + A-iBiG 1 ^ + \a~^ B^vb^G ^, (4.40) 

G\l = A-iBhG 2 Q fi +A-iBivGl% (4.41) 

G% = A-iBhG^+A-iBivG*' 2 , (4.42) 

G° 2 ° = A- 1 Bb 2 G 2 :l + 2A- 1 BvbGlll + A- 1 Bv 2 G° 2 . (4.43) 
After substituting these variables, we have effective Hamiltonian constraint without unphysical variables, 

C Q = A v \v- X - Bvb 2 - 2BbG 1 ' ! 1 . (AAA) 

So the time evolution of v is, 

= {v,C Q } = 2Bvb + 2BGl]l. (4.45) 

Then we can calculate Hubble parameter, 

Hi = £ = ^fV +2,- 1 <J] = ^[AB-\\v- 2 ] = *£p, (4.46) 

where p = 2{ J^, is again the energy density. Then we come to the same result. 

As we know the non-commutative property of operators is a key feature of quantum theory. Even though now we cast 
quantum mechanics into classical type of equations, this key feature is caught in Moyal *-product of non-symmetric 
constraints or in order symmetrization of operators in definition of symmetric constraints. 
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B. Effective Constraints of Loop Quantum Cosmology 

For k — FRW cosmology model, we consider simplified Loop Quantum Cosmology. Its Hamiltonian constraint 
operator is 

3 • ■ 1 » 2 

C = — sin(2A6)i>sin(2A6) + — — — = Ap 2 ^ 1 - BsPs, (4.47) 

47 X 2 47rG7 v 

where A = j^q^ and B — • For convenience, sin(2A6) and cos(2A6) are abbreviated as s and c respectively, and 
their expectation values are denoted as s and c. For non-symmetric constraints in effective constraint method, we 
choose the same constraint functions Cq, C v , Cb, C p . and G^ as in WDW quantum cosmology 

C Q = A V \v~ x - Bvs 2 + Av~*G%l - 2A P ^- 2 G{'° 

+A P 2 ^- 3 G 2 '° - AXBscG 1 ^ - AX 2 Bv(c 2 - S 2 )G° % 
C v = 2Ap <t> v- 1 G\l - (Ap 2 ^- 2 + Bs 2 )Gl fi - AXBuscG 1 ^ + 2ihXBvsc, 
C s = AXAp^cG ^ - (2XA P 2 v- 2 c+2XBs 2 c)G 1 1 - 8X 2 Bvsc 2 G ' 2 

-ihX{Ap\v~ 2 c + Bs 2 c), 
C m = 2Ap^- 1 Gl -(Ap 2 v- 2 + Bs 2 )G\' -4XBv S cG° i: l 

Cp = 2A P4> is- 1 G° 1 j- (Aplv- 2 + Bs 2 )Gl fi 1 -AXBvscG^\+ihAp 4> v- 1 . (4.48) 

We have checked that these constraints form a first class system. Here G^ is decoupled from other constraints. It can 
only exert some constriction on unphysical variables G 1 \, Gq'° and Go'}. Solving these constraints, we can get 

p^ = A-iBKa-^B-K^a-^ + u-^ 

—A-iB^u-hG 2 ^ + 2XA-?BicG 1 '} ) + 2X 2 A^B^vs^ic 2 - s 2 )G° ' 2 0l (4.49) 

G\'° = A-iBhG 2 '° + 2XA-^BKcG 1 1 -ihXA-iB^c, (4.50) 

G i,o = A-iBhG 1 ' 1 + 2\A-?B?iscG° '* + ^ihA-?Bh, (4.51) 

G";o = A- 1 B s 2 Gl fi +4XA- 1 BisscGl A +4X 2 A- 1 Bis 2 c 2 G°> 2 . (4.52) 

After analyzing gauge transformation generated by constraints, we can choose gauge fixing G®'® = —\ih, Gj;° = 
and Gq'\ = 0, as well as G^ = to saturate the minimum uncertainty condition G° qGq'^ — (G^'J) 2 = 

2 11 2 

Substituting relevant variables, we can express Cq in terms of G ' , G ' and G ' , obtaining effective Hamiltonian 
constraint, 

C Q = Ap\v- X - Bvs 2 - AXBscGl'l + AX 2 Bvs 2 Gvl + ^ihXBsc. (4.53) 
The imaginary term can be understood as the additional term of Moyal *-product. 

v * s 2 = ve V&&-&%) a * = vs 2 - 2ihXvsc. (4.54) 
Thus, some terms in Cq can be expressed as Moyal *-product, 

C Q = Ap\v- X -Bv*s 2 - AXBscG 1 ^ + 4X 2 Bvs 2 G°' 2 . (4.55) 
This is corresponding to the procedure of deparameterization, 

p\ = A^B^s 2 + AXuscG 1 ^ - 4AVs 2 G°;o - 2ihXvsc). (4.56) 
Then we can solve for the real Hamiltonian p^ = ±H, 

H = A-^B^{v*s + 2XcG 1 1 - 2X 2 vsG° ' 2 ). (4.57) 
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As the same way, <j> is treated as internal time, and Cq can generate time evolution. So the evolution of v is, 
v = {v, C Q } = -^(v *C q -Cq*v) = AXBvsc + 8X 2 B(c 2 - s 2 )Gj;J - 16A 3 BvscG° ' 2 - AihX 2 B(c 2 - s 2 ). (4.58) 
Now we can calculate the Hubble parameter, 

H 2 = f_ = 1^![ S 2 C 2 + AXv -l sc{c 2 _ S 2 )G M _ 8A 2 S 2 C 2 G 0,2 _ 2^-1^2 _ ^ (4 5g) 

In order to eliminate s and c, we use again the effective Hamiltonian. From Cq = 0, we can get 

s 2 = AB^p 2 ^ 2 + [-AXv^scG 1 ^ + 4A 2 s 2 Go;o + ^Xv^sc] = AB^p 2 ^ 2 + 0(h), 
c 2 = 1- AB~ 1 p 2 v- 2 -[-AXv- 1 scGl'} ) +AX 2 s 2 Gl'l + 2ihXv- 1 sc] 1 
s 2 c 2 = AB- 1 p 2 v- 2 {l-AB- 1 p 2 v- 2 ) + {l-2AB- 1 p 2 u- 2 )[-AXu- 1 scGl^ 
+4A 2 s 2 Go;o + ZihXv^sc] = Afl-y i/~ 2 (l - AB^p 2 ^ 2 ) + 0(h), 
sc = V ' AB-^p2 v -2(\ - AB- 1 p 2 v- 2 ) + 0(h), 



- 2 - 2 - l-2AB- 1 p 2 v- 2 + 0(K). (4.60) 



c 



Put them into Hubble parameter (4.59), and keep terms to order of h, we can obtain modified Friedmann equation, 



H 2 = ^^-[AB- 1 p 2 v - 2 (l-AB- 1 p 2 u- 2 )-AX 2 AB- 1 p 2 v - 2 Gl 0i 
8ttG r. If 3 ^o. 2 ^ 



-p 



where, p = 2 (2j!g^v) 2 1S tne ener gy density, and p c = 8 ^. G ^ A a is the critical energy density. 

Now we can see, when p = p c — 27r ^ 7 a G^ 2 ,, Hubble parameter reaches zero. Suppose today's cosmology can be 
described by a semi-classical state, as we evolve this state backward in the internal time, the scale of cosmology keeps 
shrinking. However we will not touch the classical singularity. Before we reach that point, we first meet the situation 
of p = p c — 27r G-y' 2 Gq,o- And at this point the Hubble parameter vanishes, the scale of cosmology reaches its minimum 
and then starts to become larger. It is corresponding to a bounce generated by quantum effects of Loop Quantum 
Cosmology. In this sense, we solve the problem of singulary, that is to say the Big Bang singulary is replaced by this 
quantum bounce. This result is different from the former one in other methods[20, 21, 25, 26] only by a term Gq'q in 
critical energy density, which is of the order of h and we can just neglect it. 



C. 3rd order effective constraints of quantum cosmology 

So far, we only consider momenta up to order of hbar, that is the terms G"'^, a + b + c + d = 2. And now we 

consider the third order, which are of order of hi, G^, a + b + c + d = 3. It emerges something new from this higer 
order system. 

In WDW cosmology, the principle constraint is (for convenience, we simplify p^ as p, and p^ as p) 



C Q = (Ap 2 ^ 1 - Bvb 2 ) 



w 



= Ap 2 v- X - Bvb 2 + Av^G^l - 2Apv- 2 G\l + Ap 2 v- 3 G 2 '° - 2BbG 1 Q } ] - BvG° ' 2 

-Av- 2 G X 2 % + 2Apv- z G 2 {l - Ap 2 v- 4 G 3 '° - BGl' 2 (4.62) 

From Cq — 0, we can get p = A^^B^vb + 0(h) So, all constraints become, 

C Q = Ap 2 v- X -Bvb 2 + A v - 1 G^l-2A^Bi v -HG\l + Bv- 1 b 2 G 2 ^l-2BbG);l 

-BvG° Q < 2 - Av- 2 Gl'° + 2AiBiv- 2 bG\l - Bv- 2 b 2 G 3 fi - BG^ 2 , (4.63) 
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C v = 2A^B2bG\' -2Bb 2 Gl'° -2BvbG 1 1 + ihBvb 



C b = 2AiBhG° 1 'l-2Bb 2 G 1 ' 1 



C p = 2A2BibG 20 



0,0 odj.2^1,0 



bG 2 {° + Bv- 


17,2,^.3,0 
^0,0 - 




- BvGl^, 


- 2BvbG^l 


- ihBb 2 






bG\l + Bv- 


WGli 




- BvG^, 


- 2BvbG {) {\ 








bG 1 ^ + Bv- 


H 2 Gl% - 


2B6Gj;J 


- BvG\% 



(4.64) 
(4.65) 
(4.66) 

C„2 = (Ap 2 i> - B\v 3 b 2 - 2ihv 2 b]) w - 2vC v - v 2 C Q 

= 2A?BhG 2 { a Q -2Bb 2 Gl° -2BisbGl] 1 , (4.67) 
C b 2 = (Aip 2 !)- 1 ^ -ihp 2 v~ 2 b\- B[vb 4 + ihb 3 ]) w -2bC b -b 2 C Q 

= 2AiBhG%-2Bb 2 Go 2 -2BisbG{\% (4.68) 
C p2 = (Ap 4 ^ 1 - Bfvb 2 ) w -2 P C p -p 2 C Q 

= 2AiBhG° 3 -° -2Bb 2 G 2 l-2BvbG° 2 % (4.69) 
C' hv = (Ap 2 b - Bu 2 b 3 ) w - bC v - vC b - bvC Q 

= 2A^B?bG\'l-2Bb 2 Gl} > -2BvbG 1 % (4.70) 
C' vb = (Ap 2 l>- B\v 2 b 3 -ihvb 2 ]) w -vC b -bC v -vbC Q 

= ihBvb 2 + 2A^B^bG\l-2Bb 2 Gll-2BvbGl'l, (4.71) 

C bv = ^{Cl + Cl b ) = UhBvb 2 + 2A^BhG\j ) -2Bb 2 Glj ) -2BvbGl% (4.72) 

C pv = (Ap 3 — B\pv b 2 — ihpvb]) w — pC v — vC p - pvCq 

= 2A^B^bG^ Q -2Bb 2 G 2 l-2BvbG\} l , (4.73) 

C pb = (A^u^b - lihfi)- 2 } - B[pvb 3 + \dhpb 2 \) w - P C b - bC p - pbC Q 

= 2A^BHG^l-2Bb 2 G\l-2BvbG%. (4.74) 

Solving for momenta, we have 

G\l = -^ihA^B^ub + A-^BhG 2 :l + A-*B*vG$, (4.75) 
G 2 >° = A-?BhG 3 ° + A-?B^G 2 'l (4.76) 



= A-^BhG^ 2 + A-^B^G % (4.77) 

Gl fi = -^ihA- l Bv 2 b + A- 1 Bb 2 G 3 a l + 2A- 1 BvbGl+ + A- l Bv 2 Gl% (4.78) 

G°;o = -^ihA- l Bvb 2 + A- 1 Bb 2 G 2 l + 2A- 1 BvbGl a + A- l Bv 2 G° % (4.79) 

G°;° - ~ihA-lBlu 2 b 2 + A-lBh 3 Gl : ° 

+3A-^B^b 2 G 2 ll + 3A-iBh 2 bG^ 2 + A~i B^ 3 G° ' 3 , (4.80) 

G{° = A-^BhG 2 ° + A-^B^Gl^-^ihA-^B^ + A-^B^G 2 % (4.81) 

= A-iBhGl'l + A-iBivG^'l + lihA-iBh + A-^B^Gli, (4.82) 

' 8 

G° 2 ° = A^B^G 2 ^ + 2A- 1 BvbG 1 a } ) + A^B^G®' 2 - h.hA^Bvb 

+2A- 1 BbG 2 1 + 2A- 1 BvG 1 ' 2 . (4.83) 
Substituting all relevant quantities into principle constraints and expressing it into a function of quantum variables, 
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finally we obtain the effective Hamiltonian constraint, 



C Q = Ap 2 ^ 1 - Bub 2 - 2BbG 1 ] 1 + '^ihBb. (4.84) 

Though we take third order momenta into acount, the resulting effective Hamiltonian acquiring no higher order 
correction. That is to say, for WDW Quantum Cosmology, there is no third order correction in effective Hamiltonian 
constraint when considering quantum effects. Also we can see that comparing with the one (4.22) in the second order 
situation, it is only different by a term —jihBb. It comes from an effect of the order of quantum operators in C^. If 
we choose C' bv , i.e. equation (4.70) as constraint instead of C^, then the effective Hamiltonian constraint is totally 
the same as the one in second order, 

C Q = Ap 2 v- X - Bvb 2 - 2B6Gj;J + ihBb. (4.85) 
But this will not affect the Hubble parameter and effective Friedmann equation, 

v = {v,C Q } = 2Bvb + 2BGl'l-jihB, (4.86) 

H ° = h = T [fe2 + 2v ~ lhG ^ i ihu ~ % (4 - 87) 

Using Cq = 0, we can derive b 2 = AB^ l p 2 v^ 2 - 2v~ 1 bG 1 ,1 + \ifrv~ x b, then 

Hl = ^[AB-Vv- 2 ]~P, (4.88) 

2 

where p = 2(iJg^v)' 1 ^ s s ^ m ^ ne ener §y density. This is exactly the same as in second order. 
Now we continue to see the situation in Loop Quantum Cosmology, its principle constraint is, 

C Q = (Atfv- 1 - Bvs 2 )w 

= Ap 2 u- 1 - Bus 2 + Au^G^l - 2Apv- 2 G\§ + V^G^o 
-AXBscGl'l - A\ 2 Bv{c 2 - .s 2 )G°; 2 - Ai/~ 2 G\l + 2A P iy- 3 G 2 {° 

-Ap 2 v- 4 G 3 '° - A\ 2 B{c 2 - s 2 )G\'l + ^\ 3 BvscG Q ' 3 . (4.89) 

From Cq = 0, we can get p = A~t-Bt-vs + 0(h). And constraints become 

C Q = Ap 2 u- 1 -Bu S 2 + Au- 1 G%-2AiB^- 1 sG 1 1 ' + Bu- 1 s 2 G 2 fi 

-AXBscGH - 4A 2 Sz/(c 2 - s 2 )G°; 2 - Av~ 2 G\l + 2A^B^v~ 2 ! sG 2 ;° 

-Bv- 2 s 2 G 3 ° - A\ 2 B(c 2 - s 2 )G^ + ^\ 3 Bv S cG° 3 , (4.90) 

C v = 2A^B^sG\ fi -2Bs 2 Gl'l - A\BvscG\ ) i + 2ih\Bvsc + Av^G 1 ^ 

-2^^B^- 1 sG 2 ;° + Bi/~ 1 s 2 G s } £ - AXBscGl'l - A\ 2 Bv(c 2 - s 2 )Gl' 2 , (4.91) 

C s = 4XAi BhcG^ - 4\Bs 2 cG 1 1 - 8\ 2 Bvsc 2 G° ' 2 - 2ih\Bs 2 c 
+2\Av- 1 cG° 2 : 1 - AXA^B^v^scG 1 ^ - 4\ 2 A^ Bh 2 G°{ 2 

+2\Bv- 1 s 2 cG 2 1 - 4\ 2 Bs(2c 2 - s 2 )G^ - %\ 3 Bvc(c 2 - 2s 2 )G°> 3 , (4.92) 
C p = 2A^B 1 isG^l-2Bs 2 G\>l-A\BvscG°{l + Av- 1 G^ 

^A^B^^sG 1 ^ + Bv- 1 s 2 G 2 1 '° - AXBscGl j, ~ 4\ 2 Bv(c 2 - s 2 )G°; 2 , (4.93) 
G„2 = (Ap 2 v - B[v 3 s 2 - Aih\0 2 sc]) w - 2vC v - v 2 C Q 

= 2A?BKsG 2 {°-2Bs 2 G 3 '°-4\BisscG 2 1 , (4.94) 
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C sv = -{C' sv + C' va ) = ih\Bvs 2 c + AXAi BhcG\ll - AXBs 2 cG 2 ' 1 - 8X 2 Bvsc 2 G^ 2 , (4.99) 



C S 2 = (Aifi)-^ 2 - 2ih\fi>- 2 sc] - B[vs 4 + 2ihXs 3 c]) w - 2sC s - s 2 C Q 

= S\ 2 AiBhc 2 G { 2 - 8\ 2 Bs 2 c 2 Gl' 2 - 16A 3 Busc 3 G° % (4.95) 
C P 2 = (Ap^- 1 - Bp 2 i)s 2 ) w - 2 P C p -p 2 C Q 

= 2A?BhGl: -2Bs 2 G 1 2 '° -4\BisscG° 2 : 1 , (4.96) 
C' sv = (Ap 2 §- Bv 2 s 3 ) w - sC v -vC a - svCq 

= 4\A^ BhcGl'l - 4:\Bs 2 cG 2 1 -8\ 2 Bi>sc 2 G 1 a ' 2 , (4.97) 
C' vs = (Afs - B[v 2 s 3 - 2ih\vs 2 c]) w - vC s - sC u - vsCq 

= 2ih\Bvs 2 c + AXA^BhcGl'l - AXBs 2 cG 2 ' 1 - 8X 2 Busc 2 G 1 ' 2 , (4.98) 

C pu — (Ap 3 — B\pv 2 s 2 — 2ih\p0sc])w — pC u — vC v — pvCq 

= 2AiBhGlj\-2Bs 2 G 2 {%-4\BisscG{ 1 , (4.100) 
C ps = (A[]5 3 j> _1 s — ih\p 3 v~ 2 c\ — B\pvs 3 + ihXps 2 c])w — pC s — sC p — psCq 

= 4\AiBhcG° 2 -l-4\Bs 2 cG{ 1 -8\ 2 Bvsc 2 G a 1 ' 2 . (4.101) 

Solved all the constraints, we have momenta 

Gil] = -h,hA-^B^vs + A-^B^sGlj ) + 2\A-^B^vcG 1 l, (4.102) 

G 2 {° = A-iBhG 3 a '° + 2XA-^BKcG 2 'l (4.103) 

G°{ 2 = A-iBhG 1 ' 2 + 2\A-iB?vcG° % (4.104) 

G^jJ = -^ih\A- 1 Bv 2 sc + A- l Bs 2 Gl^ + A\A- 1 BvscGll+A\ 2 A- 1 Bv 2 c 2 Gl% (4.105) 

G% = -jihA^Bvs 2 + A- 1 Bs 2 G 2 :l + 4\A- 1 BvscGli + 4\ 2 A- 1 Bv 2 c 2 G ' 3 , (4.106) 

G°;° = -ih\A-iBiv 2 s 2 c + A-iBh 3 G$ + 6\A-iBhs 2 cG 2 % 

+12A 2 A-?BK 2 S c 2 GH + 8\ 3 A-?Bh 3 c 3 G° < 3 , (4.107) 



G\'° = A-7B*sGl§ + 2\A-*B*i>cGl$- -ih\A-*B*vc 

+2\A-?BhG 2 'l - 2\ 2 A-?bKsG 1 2 , (4.108) 

G°{1 = A-iBhG 1 ' 1 + 2\A-?B?vcG > 2 + lihA-3Bis 

' 8 

+2XA-iBhGl' 2 - 2\ 2 A~* B?vsG° % (4.109) 
G° 2 '° = A- 1 Bs 2 Gl° + AXA^BvscGl'l + 4A 2 A- 1 Bv 2 : c 2 l G° ' 2 - ^ihXA^Busc 

+4\A- 1 BscG 2 1 + ^\ 2 A- x Bv{2c 2 - s 2 )Gj;jj - 8X 3 A~ x Bv 2 scG^%. (4.110) 
Substituting terms into principle constraint, we finally obtain the effective constraint, 

C Q = Ap 2 v~ x - Bvs 2 - A\BscGl^ + 4:X 2 Bvs 2 G° 2 + ^ihXBsc 

+4\ 2 Bs 2 G 1 „ 2 + ^\ 3 BvscG°l 3 . (4.111) 

Here we can see the above equation acquires third order momenta correction. In other words, there are higher order 
corrections in Loop Quantum Cosmology when considering quantum effects. The imaginary term here is different 
from the one in second order (4.53). This is also due to the uncertainty of the order of operators in higher order 
constraint C sv . If we use C' sv , i.e. equation (4.97), then this imaginary term comes back to the one in second order 
situation. The higher order corrections emerging from the effective Hamiltonian constraint make crutial differences 
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in Hubble parameter and effective Fricdmann equation comparing with ones in WDW Quantum Cosmology, 
v = {is,C Q } = 4\Bvsc + 8\ 2 B(c 2 - s 2 )G 1 ll-W\ 3 BisscG° ' 2 -3ih\ 2 B(c 2 - s 2 ) 



1G 



m = 



-,1,2 

3 ' 

16A 2 5 2 



IQ^BscG 1 ^ - ^-\ 4 Bv(c 2 - s 2 )G° 3 , 



9v 2 



[s 2 c 2 + AXv^scic 2 - S 2 )Gl+ - 8\ 2 s 2 c 2 G > 2 



-UhX^sc^ 2 - s 2 ) - 8\ 2 v- 1 s 2 c 2 G 1 2 - |a 3 SC (c 2 - s 2 )G° ' 3 ]. 



(4.112) 



(4.113) 



In order to eliminate s and c in above equation, we use again the effective Hamiltonian constraint. From Cq — 0, we 
have 



s 2 = AB^p 2 ^ 2 + [-AXv^scG 1 ^ + 4A 2 s 2 Go;o + ijiftAi/ -1 



,22^0,2 



SC 



+4A 2 ^~ VG^ + - A scGq'q] = AB p v~ + 0{h), 
c 2 = 1 - AB- X p 2 v- 2 - [-4Az/- 1 scGj;J + 4A 2 s 2 G°; 2 1 + ^ihXv^sc 
+A\ 2 u- 1 s 2 G];l + \\ z scG a n l\ = 1 - AB- l p 2 u- 2 + 0(h), 



2 2 
S C 



-ih\v 1 sc + 4A~i/ 's 



2,.-l„2 / ol,2 



SC = 



c 2 -s 2 



AB- x p 2 v- 2 (\ - AB- 1 p 2 v- 2 ) + (1 - 1AB- 1 p 2 v- 2 )\-A\v- 1 scG\;l + 4A 2 s 2 G^ 
3 

+ 2 

V /^ J B-V^ _2 (1 - AB- 1 p 2 v- 2 ) + 0(h) 
\-2AB- 1 p 2 v- 2 + 0(h). 



Go'fi + ^ 3 scGo 3 o] = AB- x p 2 v- 2 (\ - AB^p 2 ^ 2 ) + 0(h), 



Substitute them into Hubble parameter, we can have the modified effective Friedmann equation, 

16A 2 B 2 



9 

8ttG 



-AB- 1 p 2 v- 2 [\ - AB- X p 2 v- 2 - 4\ 2 G° 2 - iX'^G 1 ^ 



p c \ P+ 27rG^ 



{Go'fi - v x G X q q 



(4.114) 



(4.115) 



where p = 2 (2-kG~ 1 v) 2 ^ s ^ ne ener gy density, while p c = 87rG ^ A ^ is the critical energy density. The result shows that 



there is an additional momentum v Y G\ q with the order of h? . 



V. DISCUSSION 



In the above discussion, we used the method of effective constraints to obtain the effective Hamiltonian constraints 
both in WDW Quantum Cosmology and Loop Quantum Cosmology, and then used them to analyze the dynamics of 
those systems. In this method, we assumed that there exist physical semi-classical state, and based on this state we 
calculated all the expectation values and expanded them into functions of expectation values and momenta of basic 
operators. 

Since this state is general, this method should be applicable to a wild range of situations. The analyzing of the case 
of a specific type of semi-classical state, coherent state, has convinced us its reliability. Now we compare our result 
with the one from the method of path integral using a coherent state to obtain the effective Hamiltonian [27-30]. The 
resulting effective Hamiltonian in the reference [30] is, (considering the difference in denotations) 



off 



-A 



h 2 

2^ + 127rG 
2a 2 



12ttG 



-(^ 2 ) + y^ + ^ 2 + - 



(5.1) 



It contains the fluctuation 



2a' 



of p,/, corresponding to G^'q, the fluctuation ^- of v corresponding to Gq'q and the 



,-2 

fluctuation ^ 



of b which is the same as Gq'q. Since they choose the coherent state as the semi-classical state, 
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the corelation function between 6 and v is itself vanishing. Thus there is no such term appearing in the effective 
Hamiltonian. In this effective Hamiltonian constraint it contians a fluctuation term of p^, while in our method of 
effective constraints, we have solved the corresponding term G" q and put it in terms of other functions. If we keep 
this term, then the effective Hamiltonian constraint becoms, 

C' Q = Aplv- 1 - Bvb 2 + Au^G ^ - Bv-WgI'^, - 4B6Gj;J - BvG ^ + ihBb. (5.2) 

Then we can obtain 

J^ s = -p\ - G° 2 ° + 127rG(V 6 2 + 6 2 G 2 ; I + 4i/6G5;J + v 2 G% - ihvb) 

= ~v\ - G% + \2itgL * (vb 2 ) + 6 2 G 2 ; I + 4i/&G$;J + ^ G ofi) ■ ( 5 - 3 ) 



From above equations, we can see these two results match to each other without higher order terms. 

For comparison, we change the variables both in WDW Quantum Cosmology and Loop Quantum Cosmology, 

2|p|* 27rG 7 VA 

where they satisfy the basic Poisson bracket {b, v} = 1. In Loop Quantum Cosmology, the constraint operator is, 

3 • ■ 1 v 2 

C = sin(26)z>sin(26) + — — ^- = Ap 2 ^ 1 - Bsus, (5.5) 

47 47rG7 v 

where A — , B = In the paper [30], they use the coherence state to get the Hamiltonian, 

4ft 2 , 2h 2 \ 4ft 2 



^off = -4-^2 + 37rG 



(v 2 + d -) (sin 2 (26)(l - ^-) + ±L) - ifo/sin(46)(l 

= -Pi £5 +**G[»* (-sin 2 (26)(l - if)) + ^v 2 + t sin 2 (26)(l - + ft 2 

(5.6) 

where it contains the fluctuation ^ of p<j, corresponding to G^'g, the Auction ^- of v corresponding to Gq'q and the 
fluctuation 2 Jj- of b which is the same as Gq'q. Additionally it includes a higher order correction term ^5- of the 
variable sin 2 (26). We suppose that this term is ignorable, then 

,2 d 2 



iTcff = -v\ - ^2 + 37rG 



(i/ 2 + y) (sin 2 (26) + ^L) - zftz,sin(46) 



2 



d 2 . 2ft 2 



^* (^sin 2 (26)) + — sin 2 (26) + — v 2 + h 2 . (5.7) 
2 a 2 

Because they use the coherence state as the semi-classical state, the corelation function between 6 and v is itself 

2 

vanishing. In the method we introduced, keeping G ; , the hamiltonian constraint is, 



C' Q = Ap 2 v~ x - Bvs 2 + Av~ 1 G°2fi - Bv~ 1 s 2 G 2 ; () 

-SBscGlll - ABv{c 2 - s 2 )G a 2 + 2ihBsc. (5.8) 

Then we can obtain the hamiltonian, 

Jt^ s = -p 2 -G a 2 : + 37TG^iy 2 s 2 + s 2 G 2 a + 8iyscGl:} ) + A v 2 (c 2 - s 2 )G ' 2 - 2ifii/sc) 

= -P 2 - G";o + 3ttG(^ * {vs 2 ) + s 2 G 2 ° + 8i/scGj;J + 4^ 2 (c 2 - s 2 )G° fy . (5.9) 

So ignoring the higher order terms, there is only one term v 2 {c 2 — s 2 )Gq'q which is different to the one v 22 ^r in our 

method. This difference comes from the mismatching between 2 % r in (5.7) which is the fluctuation of sin 26 and Gq'q 
in (5.8) which is the fluctuation of 6. The reason underlining this diffcrccc still needs further investigation. 

In conclusion, it shows that the method of effective constraints we used to handle quantum constraint systems, 
solving unphysical variables and keeping only physical variables to obtain the effective Hamiltonian constraint is 
reasonable. Additionally, the property that the classical big bang singularity is replaced by the quantum bounce in 
Loop Quantum Cosmology is in fact reliable. 



17 

Acknowledgments 

This work is supported by NSFC (No. 10975017) and the Fundamental Research Funds for the Central Universities. 



[1] A. Ashtekar and T. A. Schilling, Geometrical Formulation of Quantum Mechanics, CGPG 97/6-1, [math-ph/0511043] 
[2] M. Bojowald and A. Skirzewski, Effective Equations of Motion for Quantum Systems, Rev. Math. Phys. 18 (2006): 713-745, 
[math-ph/0511043] 

[3] M. Bojowald and A. Skirzewski, Quantum Gravity and Higher Curvature Actions, Int. J. Geom. Meth. Mod. Phys. 4 

(2007) : 25-52, [hep-th/0606232] 

[4] M. Bojowald, B. Sandhofer, A. Skirzewski and A. Tsobanjan, Effective constraints for quantum systems, Rev. Math. Phys. 

21 (2009): 111-154, [arXiv:0804.3365] 
[5] H. Feshbach and F. Villars, Elementary Relativistic Wave Mechanics of Spin and Spin 1/2 Particles, Rev. Mod. Phys. 

30 (1958): 24-45 

[6] J. B. Hartle and D. Marolf, Comparing Formulations of Generalized Quantum Mechanics for Reparametrization-Invariant 

Systems, Phys. Rev. D 56 (1997): 6247-6257, [gr-qc/9703021] 
[7] M. Bojowald, Large scale effective theory for cosmological bounces, Phys. Rev. D 75 (2007): 081301(R), [gr-qc/0608100] 
[8] M. Bojowald, Dynamical coherent states and physical solutions of quantum cosmological bounces, Phys. Rev. D 75 (2007): 

123512, [gr-qc/0703144] 

[9] M. Bojowald, H. Hernandez, and A. Skirzewski, Effective equations for isotropic quantum cosmology including matter, 
Phys. Rev. D 76 (2007): 063511, [arXiv:0706.1057] 
[10] M. Bojowald, Quantum nature of cosmological bounces, Gen. Rel. Grav. 40 (2008): 2659-2683, [arXiv:0801.4001] 
[11] M. Bojowald, Harmonic cosmology: How much can we know about a universe before the big bang?, Proc. Roy. Soc. A 464 

(2008) : 21352150, [arXiv:0710.4919] 

[12] M. Bojowald and R. Tavakol, Recollapsing quantum cosmologies and the question of entropy, Phys. Rev. D 78 (2008): 
023515, [arXiv:0803.4484] 

[13] M. Bojowald and T. Strobl, Poisson Geometry in Constrained Systems, Rev. Math. Phys. 15 (2003): 663-703, [hep- 
th/0112074] 

[14] A.Komar, Constraints, Hermiticity, and Correspondence, Phys. Rev. D 19 (1979): 2908-2912 
[15] CRovelli, Quantum Gravity, Cambridge: Cambridge University Press, 2004 

[16] A. Ashtekar, J.Lewandowski, Background Independent Quantum Gravity: A Status Report, Class. Quant. Grav. 21 (2004): 
R53R152, [gr-qc/0404018] 

[17] T. Thiemann, Modern Canonical Quantum General Relativity, Cambridge: Cambridge University Press, 2007 
[18] M. Han, Y. Ma, W. Huang, Fundamental Structure of Loop Quantum Gravity, Int. J. Mod. Phys. D 16 (2007): 13971474, 
[gr-qc/0509064] 

[19] A. Ashtekar, M. Bojowald, J. Lewandowski, Mathematical Structure of Loop Quantum Cosmology Adv. Theor. Math. 

Phys. 7 (2003): 233268, [gr-qc/0304074] 
[20] A. Ashtekar, T. Pawlowski, P. Singh, Quantum Nature of the Big Bang: An Analytical and Numerical Investigation, Phys. 

Rev. D 73 (2006): 124038 [gr-qc/0604013] 
[21] A. Ashtekar, T. Pawlowski, P. Singh, Quantum Nature of the Big Bang: Improved Dynamics, Phys. Rev. D 74 (2006): 

084003. [gr-qc/0607039] 

[22] A. Ashtekar, A. Corichi, P. Singh, Robustness of Key Features of Loop Quantum Cosmology, Phys. Rev. D 77 (2008): 
024046, [arXiv:0710.3565] 

[23] A. Ashtekar, P. Singh, Loop Quantum Cosmology: A Status Report, Class. Quant. Grav. 28 (2011): 213001, 
[arXiv: 1108.0893] 

[24] J. Yang, Y. Ding, Y. Ma, Alternative Quantization of the Hamiltonian in Loop Quantum Cosmology II: Including the 

Lorentz Term, Phys. Lett. B 682 (2009): 1-7, [arXiv:0904.4379] 
[25] V. Taveras, Corrections to the Friedmann Equations from Loop Quantum Gravity for a Universe with a Free Scalar Field, 

Phys. Rev. D 78 (2008): 064072, [arXiv:0807.3325] 
[26] Y. Ding, Y. Ma, J. Yang, Effective Scenario of Loop Quantum Cosmology, Phys. Rev. Lett. 102 (2009): 051301, 

[arXiv:0808.0990] 

[27] A. Ashtekar, M. Campiglia, A. Henderson, Path Integrals and the WKB Approximation in Loop Quantum Cosmology, 

Phys. Rev. D 82 (2010): 124043. [arXiv: 101 1.1024] 
[28] L. Qin, H. Huang, Y. Ma, Path Integral and Effective Hamiltonian in Loop Quantum Cosmology, Submitted to Gel. Rel. 

Grav. [arXiv:1102.4755] 

[29] L. Qin, G. Deng, Y. Ma, Path Integrals and Alternative Effective Dynamics in Loop Quantum Cosmology, Commun. 

Theor. Phys., 57 (2012): 326. [arXiv: 1102.4755] 
[30] L. Qin, Y. Ma, Coherent State Functional Integrals in Quantum Cosmology, Phys. Rev. D 85 (2012): 063515. 

[arXiv: 11 10.5480] 



